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Abstract — In this paper we investigate a multi-source LDPC 
scheme for a Gaussian relay system, where M sources com- 
municate with the destination under the help of a single relay 
(M — 1 — 1 system). Since various distributed LDPC schemes 
in the cooperative single-source system, e.g. bilayer LDPC Q) 
and bilayer multi-edge type LDPC (BMET-LDPC) |2|, have been 
designed to approach the Shannon limit, these schemes can be 
applied to the M—l — 1 system by the relay serving each source 
in a round-robin fashion. However, such a direct application is 
not optimal due to the lack of potential joint processing gain. 
In this paper, we propose a network coded multi-edge type 
LDPC (NCMET-LDPC) scheme for the multi-source scenario. 
Through an EXIT analysis, we conclude that the NCMET-LDPC 
scheme achieves higher extrinsic mutual information, relative to 
a separate application of BMET-LDPC to each source. Our new 
NCMET-LDPC scheme thus achieves a higher threshold relative 
to existing schemes. 

Index Terms — Multi-source LDPC, network coding, network 
capacity, extrinsic mutual information. 

I. Introduction 

Low-density parity-check (LDPC) codes have been shown 
to approach theoretical capacity limits for single link com- 
munication channels [1|. Recently, distributed LDPC for co- 
operative communications has attracted much attention. The 
work of 12 first explored the the use of bilayer LDPC codes 
within the cooperative single source channel (1 — 1 — 1 system), 
where full-duplexing relay is used. Although bilayer LDPC is 
carefully designed to approach the system capacity [4], the 
performance is decreased as the capacity gap between the 
source-to-relay channel and the source-to-destination channel 
becomes larger. In 0, multi-edge type LDPC code has been 
utilized to address this problem 0. The works of |7|, |8| 
consider more practical issues in the 1 — 1 — 1 system, such as 
the use of Rayleigh fading channels and half-duplexing relays. 

However, the above studies on distributed LDPC codes are 
all limited to the triangle model, which contains only one 
source. In this paper, we investigate a network coding |9| 
based LDPC codes designed for the cooperative uplink system 
with multi-source and one relay (M — 1 — 1 system) as 
shown in Fig. Q] Based on existing methods of distributed 
LDPC design in the triangle model, an intuitive thought is 
that the relay serves the sources in a round-robin fashion, 
optimizing the distributed LDPC for a single source in each 
round. Unfortunately, such a direct application is not optimal 
for the following reasons. 




Fig. 1. Multi-source relay system. 

1) The check digits produced by the relay in the m-th, 
(m = 1, • • ■ , M), round are only based on the codeword 
from the source s m , which is highly correlated with the 
original check codes already produced by s m . 

2) Network coded check digits enable the joint decoding of 
all sources' data and thus each source can obtain more 
extrinsic mutual information from the other sources. 

3) The code profile optimization executed by the relay aims 
only at approaching the capacity for the single source, 
which leads to a sub-optimal outcome for a multi-source 
system in terms of network capacity. 

We note that iflOl has studied the network coded LDPC 
in a multi-source system with fading channels. But the code 
design in 1 1 1 is not optimal. [ 1 1 1 also proposes a joint bilayer 
LDPC scheme in the M—l—1 system. However, this scheme 
is constrained by the drawbacks of bilayer LDPC. So it cannot 
deal with the problem where the capacity gap between source- 
to-relay channel and source-to-destination is large. Also, ifTTl 
only considers a scenario where all channel capacities in the 
model are equal. In this paper, we propose a network coding 
based multi-edge type LDPC, which we refer to as NCMET- 
LDPC, for the M — l—1 model which addresses these issues. 
In our analysis of NCMET-LDPC we utilize extrinsic mutual 
information (EMI) transfer, rather than density evolution (DE) 
used in IfTTl . which aids in understanding the code design 
issues better. Through the EXIT chart, we see that the network 
coded parity check digits produced by relay provide more EMI 
for each source, relative to a separate round-robbin scheme. 

II. System Model and Preliminaries 

Consider a Gaussian relay system with A/-sources, 1- 
relay and 1 destination as shown in Fig. [T] where sources 
si, ■ ■ ■ ,sm ttansmit information to the destination d simul- 
taneously with the help of a full-duplex relay r. We assume 



2 



| Xlm~MO,Plm) Z2m~N(0,Nlm+N2m) 
1 



Y\m '. X\m 




Xm~N(0,Pm) 



Fig. 2. The Gaussian relay system working at frequency band / m : X m is 
power constrained to P m ; Xi m is power constrained to P\ m - 
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all sources are randomly distributed around the relay. Suppose 
that the m-th source, s m , transmits its information in the fre- 
quency band f m , and r can receive and transmit at fx, ■ ■ ■ , /a/. 
All the frequency bands are assumed to be orthogonal. With 
these constraints the multi-source system can be viewed as 
M independent parallel 1 — 1 — 1 systems. Within each of 
these channels, a bilayer LDPC scheme can be utilized to 
approach the 1 — 1 — 1 system capacity J4J. 

Let us briefly review the achievable rate and code design 
in the 1 — 1 — 1 system (2]. Without loss of generality, 
we focus firstly on a specific 1 — 1 — 1 system. This is 
formed by selecting a specific value of m, and constructing 
a model composed of s m , r m and d, where r m is the part of 
r operating at f m . In Fig. [2] X m is the signal transmitted by 
s m , which has the average power P m , and X\ m is the signal 
transmitted by r m , which has the average power P\ m . The 
binning scheme is used to achieve the capacity in a Gaussian 
degraded relay channel [4 |. In this binning scheme, s m divides 
its total power P m into a fraction aP m for the new codeword 
u>i, and a fraction (1 — a)P m for the bin index fa of the 
previous codeword u>i-x- So in the i-th time slot, X m is the 
superposition of u>i and fa, which will be received by both r m 
and d. Since r m has successfully decoded LUi-i in the previous 
time slot, u) i will be successfully decoded at the relay with a 
rate no more than 

aP„ 



(1) 



Meanwhile, r m is transmitting Xx, m (fa) to d. Thus, d receives 
the interferential signal composed by X m and Xxm in the 
frequency band f m . By successive interference cancellation, 
d firstly treats X m as noise so as to extract the bin index fa 
with a rate no more than 



{VPT^+V(l-a)P m ) 2 
aP m + N lm + N 2m 



(2) 



Then combining with fa, the decoding of uJi-i at d will be 
successful with a rate no more than 



R m = - log ( 1 + — 

- 2 H N lm + N 2m 



(3) 



Combining the above three equations, we get the overall rate 
for the Gaussian relay channel working at f m as 



R m = maxmin{i?™ R? + R m }. 



(4) 



So to maximize R m , we let i?™ = R™ + R™ by adapting 
power allocation a. 

For the practical code design, we can utilize the bilayer 
LDPC scheme for each s m (m = 1, ■ • ■ , M) to approach i?" 1 



Fig. 3. Comparison of the two LDPC code schemes, (a) The conventional 
bilayer LDPC code design for each source individually; (b) The proposed 
multi-edge type LDPC with joint processing. 



and R™ |2]|. For example, the code structure is divided into the 
lower and the upper graphs and consists of kxm an d check 
nodes, respectively. Both type of check nodes are connected 
to the same n m variable nodes. According to [2], in the design 
of such codes, one should firstly determine the optimal LDPC 
code corresponding to the lower graph to achieve rate i?™, 
and then search over the whole bilayer graph in order to 
find a good LDPC code that approaches the rate R™. So in 
the multi-source case, a straightforward technique of practical 
code design is to perform bilayer LDPC individually to each 
source. As previously mentioned, such individual processing 
does not benefit from any joint processing gain at the relay. In 
the following section, we focus on the network coded multi- 
edge type LDPC code for the multi-source system. 

III. Multi-edge Type LDPC for Multi-source 
System 

A. Multi-edge Type LDPC Codes 

The principle of multi-edge type LDPC is to introduce more 
than one edge type to the Tanner graph J5), where the graph 
ensemble is specified through two polynomials, one associated 
to variable nodes and the other associated to constraint nodes. 
The two polynomials are given by 

v(r,x) = ^w b ,dr b x d and /i(x) = ^MdX d , (5) 

where d = (dx, ■ ■ ■ , d ni ) is a multi-edge degree and x = 
(xx, ■ ■ ■ , Xm) denotes variables. Similarly, b = (bo, ■ ■ ■ , b nT ) 
is a received degree, and r = (rrj , • • ■ , r„ T ) denotes variables 
corresponding to received distributions. We use x d to denote 
ll™=i x i an d rb to denote n"=o ^ ■ F° r more details about 
multi-edge type LDPC codes one can refer to 0. 

B. Design of Multi-edge Type LDPC 

In this subsection, we propose a novel multi-edge type 
LDPC scheme in the M — 1 — 1 system. In the new coding 
scheme, we do not change the power allocation of each 
transmitter (i.e network capacity is not altered), but propose a 
new coding scheme to approach the network capacity. Fig. [3] 
compares the proposed scheme with the conventional bilayer 
LDPC code. In the conventional scheme, the additional k2m 
check digits at r m are only produced from s m 's frame and 
transmitted at f m . However, in the proposed scheme, the extra 
fe m bits of information are co-produced from all the sources' 
frames and randomly distributed in all the frequency bands. 



3 




km "' km 



Fig. 4. Tanner graph and multi-edge type in proposed scheme. 



r jointly processes all sources' information and produces 
&2 = X)m=i^2m parity check digits, which can be seen as a 
super parity check block. The Tanner graph of the multi-edge 
type LDPC is shown as Fig. @] where we represent the edge 
types as £ with different subscripts. The edge of the lower 
graph of s m is denoted as £± m and the edge of the upper 
graph is denoted as E^m- We also assume the frames from all 
the sources have the same length, i.e., m = • • • = njj = n. 

The multi-edge type LDPC code design for an M — 1 — 
1 system begins with optimizing the lower Tanner graph in 
Fig. @]at rate R\ , • • • , R¥ for s\ , ■ ■ ■ , sm, respectively, which 
follows the conventional methodology of single link LDPC 
codes. So the lower graph ensemble of the multi-edge type 
LDPC codes for s m is represented by 

d„i, m 

v m (r,x)=ri v [om<ii, m ] x t™ > 



(6) 



/i, 



di m = l 



where [0, 1] is the vector b, and [di >m ] is the vector d of ([5]). 
For vector b, all variable nodes in the codeword are transmitted 
through the source-to-relay channel (b\ = 1) at rate i?" 1 and 
there are no punctured variables (bo = 0) in the codeword. 
Vector d contains only one element since there is only one 
edge type. d\ >m is denoted as the degree of variable nodes 
and check nodes with the maximum value d v i m and d c i, m , 
respectively. The quantity Ub.d^ is the the number of variable 
nodes of type (b, d) and /idn is the number of check nodes 
of the type d in the graph. The code rate of the lower graph 
ensemble for s m is 



E 



M[d 1 



(7) 



So the sum of code rate for the lower graph ensembles for all 
the sources is 



M 



M dcl.m 



In the next step, we design the overall graph considering the 
relay jointly processing all of the sources' information. Relay r 



will transmit additional &2 = X)m=i^ 2 »™ P^ity check digits to 



d within the relay-to-destination channel capacity X)m=i^T' 
in which R" l n bits are allocated to s m . So the variable and 



check nodes' polynomials for all sources s m in the overall 
graph can be written as 

d v l,m d v2 ,m 

V m (r,X)=n Y E V {0A[di, m 42, m ] X l,m x 2,m ( 9 ) 
di m =ld 2 m =0 



del. 



i( x ) = E M[di,m,0,-,0]*l, l m + E ^[0.^2,1. 



,d 2 ,M] X 2,m 



dl m=l 



These relations mean that for the overall graph ensemble, all 
variable nodes in the codeword are transmitted through the 
source-to-destination channel (b\ = 1) at rate R™, and there 
are no punctured variables (po = 0). Vector [di >m) <^2,m] m ® 
represents the number of the sockets of the two edge types £ i OT 
and £% m in a variable node with c?i. m and e?2,m respectively. 
Vector [di. m ,0, ••• ,0] in (0 represents the number of the 
sockets of the edge type £ \ m in a check node of the lower 
graph with c?i. m . Since the check nodes in the lower graph 
are only connected to £\ m , the number of the sockets of other 
edge types is all zero. Vector [0, ^2,1, • • • , da m] m me second 
equation of (O represents the number of the sockets of the 
edge types £21, • • • , £im m a check node of the upper graph 
with da i, • • • , d2,M- Since the check nodes in the upper graph 
are not connected to £\ m , the number of the sockets of this 
edge type is zero. Note that for s m we have 



d c2 ,„ 



R m <l- Y M[d lim ,0,-,0] - E E ^[O.da.l.- 
di TTi— 1 d 2 m — 1 ~d2 jTi 



where 



d c2 



dc 2 



E - E ■•• E E ••■ E 



(10) 



(11) 



1=0 



1=0 d 2 , m + i=0 



f=0 



Since the k% parity check digits are shared by M types of 
edges, their contribution to s m is 



d c2 ,„ 



R? 
Then we get 



E E 

d 2 .m = l ~d 2 .n 



M[0,d 2 ,i,--- ,d 2 ,M]"2,m 



(12) 



dd, 



R"^ — 1 ~ E ^l.m.0,-,0] _ -ST 



(13) 



di m = l 



So the code rate of the overall graph ensemble can be 
computed as follows. 



r- = m-J2\ E /*[*.«,<>,••• ,o]+izr 



(14) 



. di m =l 



We optimize the whole system by regarding all the sources' 
frames as a super block, which accesses the relay with code 
rate i? + for the lower graph, and accesses the destination with 
code rate i?_ for the overall graph. Besides (0, ( [12] ) and ( Tl~3b . 
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there are several constraints that should be satisfied for each 
s m as follows; 

d„ 2>m 

U [0,l],[d ltm ] = E V V>A],[<ll,m,d2, m h 

d,2,m=a 

M[di, m ] = M[di, m ,0,-" ,0] 
d c i m d v i m 

E 



dl. 777 = 1 



(15) 



E 



1^2.777. — 



E M[0,d 3 ,i,-,<l3,Af] a 2,m 

d 2> m = l "^2,777 

d„i, 



E E W [0,l],[rfi. m ,ci2, m ] d 2,m- 



ll,m = l <l2, m = 



To deduce the average extrinsic mutual information of each 
edge type, we characterize the code ensemble of s m by the 



degree distribution A 



[di, m ,d 2im ] 



for i = 1,2, 



A 



u [0,l],[<ii, m ,<J2,m]"i-- 



Z^d 1 , ! = 1 Z^d 2 , i= o t; [0,l],[di, ! ,d2.!]"l,i 



(16) 



This defines the percentage of £j m type edges connected to the 
variable nodes with di >m edges in £\ m , and d2, m edges in £2™- 
We also define another two types of degree distribution. One 
is ph™ p which denotes the percentage of £\ m type edges 
connected to the check nodes in the lower graph with di iTn 
edges, i.e., 

P[i lirr i = — . (17) 



and the other is p^J™ 2 1 dj i, which denotes the percentage 
of £ 2 m type edges connected to the check nodes with the edge 
vector [0, d 2 ,i, • • • , c^.m], i.e., 

2,m _ / i [0,d 2 ,i,"' ,d 2 ,M] c '2,m 

^[0,d 2 ,i,--- ,<J 2 ,m] 



E«C2,771 V-V 
d 2 i=o2j~d 2 ^[O.di, 



(18) 



'2./ 



IV. Performance Analysis and Code Design 

In 0, Density Evolution (DE) is applied to the bilayer 
LDPC code profile optimization. Due to the fixed degree 
of the check nodes in both the lower graph and the upper 
graphs, the complexity of DE is tolerable. However, the fixed 
degree deteriorates the system performance as mentioned in 
the introduction. So in the optimization of the multi-edge 
type LDPC code, we exploit the extrinsic information transfer 
(EXIT) functions 11121 . fl3l . |[T4| to reduce the code searching 
complexity. This will likely lead to better code profiles by 
canceling the constraints on the check nodes degree as in 
bilayer LDPC code. 

We denote the variable nodes set associated with the code- 
word bits of s m as V m , and the check nodes set associated with 
the parity check digits of s m from the lower graph as Ci, m . 
The shared check nodes set in the upper graph is denoted as 
C*2- Since V m is connect to two edge types, i.e., C\_ m and 
C2, there are four types of mutual information (MI) defined 
as follows Ifl4l . 



/E„(l,m): The MI between the message sent from V m to 
Ci m and the associated codeword bit, on each edge in the 
edge type £ lm connecting V m to C\^ m . 

Iev(2i m )' The MI between the message sent from V m to 
C2 and the associated codeword bit, on each edge in the edge 
type £ 2m connecting V m to C 2 . 

Je c (1, m): The MI between the message sent from C\ >m to 
V m and the associated codeword bit, on each edge in the edge 
type £ lm connecting Cx >m to V m . 

7b c (2, m): The MI between the message sent from C2 to 
V m and the associated codeword bit, on each edge in the edge 
type £ 2 m connecting C 2 to V m . 

Note that the extrinsic MI on an edge connecting V m 
to Ci tm (oT C2), at the output of the variable node, 
is the a-priori MI for Ci. m (or C 2 ), i.e., /^(Lm) = 
1,40(1, m) (or Iev(2> fn) = Iac{^, m ))- Similarly, the extrinsic 
MI on an edge connecting Ci )Tn (or C2) to V m , at the 
output of the check node, is the a-priori MI for V m , i.e., 
I Bc (l,m) = I Av (l,m) (or I Ec (2,m) = I Av (2,m)). Then 
we have the iterative process as follows. 

1. Variable nodes to check nodes update. The mean of 
the extrinsic MI on an edge type £\ m connecting V m to Ci >m , 
at the output of the variable node in the l-th iteration is 



d v l,m d v 

(!."»)= E E ((di,m-l) J-\l A l l(l, m )) 

^1,777= 1 ^2,777=0 
2 



l2,m 



J 



- 1 (/W(2,m))] +[j-\l ch (m))] 2 U^ m ^ m] 



(19) 



Also, the mean of the extrinsic MI on an edge of £2™ 
connecting V m to C 2 , at the output of the variable node in 
the l-th iteration is 



rf 2 ,m=ldl,m=l 



(2."»)= E E ( rf 2,7,7-l) J- 1 (41(2, m)) 



J 



-\lt(l,m))]\[j-\l ch (m))f) \\£ m4 . 



(20) 



We can now get the MI in the l-th iteration as I l Ev (l,m) = 
J te'ftm)) and I l Ev (2, m) = J U<$%m)\ . 

2. Check nodes to variable nodes update. The update 
from check node to variable nodes is more complicated. We 
give the approximation according to 1T31 . The extrinsic MI on 
an edge type £\ m connecting C\^ m to V m , at the output of the 
check node in the l-th iteration is 

7«(l,m) = l- 



\ dl,m = l 



(d llTn -l) J- x {l-I { Z{l,m)) 



l.m 

^1,777 



(21) 



The extrinsic MI on an edge type £2771 connecting C2 to 
V m at the output of the check node in the l-th iteration is 
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more complicated as more than one source participates in the 
generation of Ci- We have 

/£(2,m) = l- 



J 



^2. 777 = 1 ~d 2 



(d2,m-l) J- 1 (l-4o( 2 » m )) 



M 

£*.« 

m' — 1 



J-i(l-7«(2,m')) 



2,m 
^[0,<*2,i, 



/ 

(22) 



In each iteration process, we make lAv{h m) — lEc(i> m ) and 
lA c (h m ) — lEv(i>™) f° r « = 1,2. Ichiro) is determined ac- 
cording to the SNR of the m-th source-to-destination channel. 
At the end of the iteration, the MI between the V m and the 
associated codeword is 



J(m) = J 



"ill, 77! "7j2,777 

£ £ 

^l,m=1^2,m=0 



J-i(I^(l,m)) 



2 

l2 \ \ . 



J-i(l£(2,m))j +[^- 1 (i-ch(m))r)2JA{3 1 , m|(faira] 

i=l 

(23) 



where 



U[0,l],[di, m ,d 2 , m ] 



(di, m + rf 2 



(24) 

According to the design of bilayer LDPC code, we first fix 
the lower graph codes of all sources, choosing the optimal 
point-to-point LDPC code to approach capacity Then 
we optimize the overall graph to approach the capacity of 
the whole system. The code optimization involves finding M 
variable node degree distributions, i.e., W[o,i],[di m ,d 2 m ] ^ or Sm 
and a check node degree distribution fj,m,d 2 i,--- ,d 2 M \- The 
optimization problem is concluded as minimizing the SNR 
of the whole system, which is a dual problem of maximizing 
the system threshold. We can therefore write 



1 



maximize a 



.sy.s 



V t^TF i<? M ) 2 ( 25 ) 

subject to I(m) — > 1, for m = 1, • • • , M. 

V. Numerical Results 

We choose the 2-sources case to illustrate the design. The 
first source has the following capacities: R\_ = 0.7, B}_ = 0.5 
and R\ = 0.2. The second source has the following capacities: 
R\ = 0.58, R 2 _ = 0.38 and R\ = 0.2. Applying the proposed 
method, we get the code profile shown in Table Q] Note that 
o\_ and al_ in Table U are the thresholds for R\ and R\, 
respectively. Also note that the threshold deduced by EXIT is 
always larger than that deduced by DE, which has been proved 
by lfT6l . So a\_ and a 2 _ in the Table U are larger than the exact 
values, which are searched by criterion ( f25l l. 



Source 1 



Source 2 



Variable Node Distribution 



f[0,l][dl,l.*2,ll 


£n 


£21 


"[0,1] [dl,2,d 2 ,2] 


£12 


£22 


0.244225 


2 





0.3289203 


2 





0.1531552 


2 


1 


0.0772109 


2 


1 


0.0209422 


2 


7 


0.0531292 


2 


2 


0.1930021 


3 





0.145309 


3 





0.138759 


3 


3 


0.0149215 


3 


1 








0.123802 


3 


2 








0.0286741 


3 


14 


0.058304 


6 





0.0346943 


6 





0.0109062 


6 


7 


0.0101216 


6 


2 


0.000131148 


6 


21 








0.05680712 


7 





0.092595 


7 





0.047728 


7 


2 


0.0297043 


7 


7 


0.0556525 


20 


3 


0.0165257 


20 





0.0203505 


20 


7 


0.00437642 


20 


1 








0.0400163 


20 


3 


Check Node Distribution in Lower Graph 


^[di.l] 




^[dl. 2 ] 


£12 


0.3 


15 


0.42 


10 


Check Node Distribution in Upper Graph 


M[0,d 2 ,i,d 2i2 ] X d 2 i + d 2 2 


£21 


2 

M[o,d 2il ,d 2>2 ] x ch ± + d2 2 


£22 


0.4 x 0.5 


3 


0.4 x 0.5 


3 


Code Rate and Threshold 


R\ = 0.7, R L _ = 0.5 


R\ = 0.58, R± = 0.38 


a\ = 0.722955, a L _ = 0.970555 


<t\ = 0.859273, ai = 1.189900 



TABLE I 

The node distribution of the NCMET-LDPC in two-source case. 



First of all, we determine the distributions of £n and 
£12 in the lower graphs for si and s 2 to approach the rate 
R\. and R\, respectively. They are designed as the single 
link LDPC codes and can be directly obtained from |15|. 
Our main task is to find out the optimal distribution of £21 
and £ 22 , and the corresponding U[o,i][c2 lil ,d 2 ,i]> w [o,i][di, 2 ,d 2 , 2 ] 
and ^[o,d 2 !,d 2 2 ]- By adopting the searching criterion of (|25T > 
and EXIT curves fitting, we get the three elements with 
the thresholds a x _ = 0.970555 and <r 2 _ = 1.1899 for s x 
and S2, respectively. The variable node distributions of the 
two sources, U[o,i][d!i,i,da,i]> w [o,X][di,3,d2,2]> and corresponding 
degrees are shown in Table U The check node distribution in 
the upper graph at the relay is jLtm^j m ,d 2 m ]> which has only 
one distribution as m .d 2 m ] = R- + R- = 0.4. Each 

check node of the relay has the degree 6, half of which are 
allocated to s% and the other half are allocated to s%. This is 
reasonable since i?L = R 2 _ = 0.2. 

Fig. |5] shows the EXIT charts for the 4 edge types at the 
system threshold. Fig. [6] and Fig. [7] show the EXIT chart at 
a\_ = 0.969555, a 2 _ = 1.1909 and a x _ = 0.971555, a 2 _ = 
1.1889, respectively. Note that in the three figures, we adopt 
the same value of the summation of a\_ and a 2 _ . However, the 
system threshold in Fig. [6] is larger than that of Fig. [5] and 
the system threshold in Fig. [7] is smaller than that of Fig. [5] 
So we can see that the EXIT curves in Fig. [6] intersect at a 
value smaller than 1, which means that the iterative decoding 
at destination will not converge eventually. The EXIT curves 
in Fig. Q as well as Fig. [5] intersect at 1, which means the 
decoding at destination will succeed. Finally, we investigate 
the separate LDPC code without network coding. In this case, 
each source is only connected to half of the check nodes at 
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Fig. 6. EXIT curves with ai = 0.969555 and o 2 _ = 1.1909. 
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Fig. 7. EXIT curves with <y x _ = 0.971555 and o 2 _ = 1.1889. 



the relay, and each check node at relay only contains one 
edge type, either £2,1 or £2.2- On the other hand, we keep the 
variable node distributions of each source unchanged. Fig. [8] 
shows the EXIT curves of the two edge types connected to 
s\. Obviously, both edge types cannot converge at 1. So we 
conclude from this case that the separate LDPC code obtains 
less extrinsic mutual information from the check nodes at relay 
than the proposed coding scheme. 

VI. Conclusion 

In this paper, we investigate a network coded LDPC design 
in the multi-source scenario. We apply the multi-edge LDPC 
to the system and execute the EXIT analysis. We conclude 
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Fig. 8. EXIT curves fit for two sources at the system threshold a\_ 
0.970555 with rate R 1 = 0.5 and cr 2 = 1.1899 with rate R 2 = 0.38. 



that each source achieves more extrinsic mutual information 
due to joint processing at the relay. Therefore, our scheme 
delivers better performance compared to traditional schemes 
that do not utilize joint processing at the relay. 



References 



[i] 



T. J. Richardson, M. A. Shokrollahi, and R. L. Urbanke, "Design of 
capacity-approaching irregular low-density parity-check codes," IEEE 
Trans. Inf. Theory, vol. 47, no. 2, pp. 619-637, Feb. 2001. 
P. Razaghi and W. Yu, "Bilayer low-density parity-check codes for 
decode-and-forward in relay channels," IEEE Trans. Inf. Theory, vol. 53, 
no. 10, pp. 3723-3739, Oct. 2007. 
[3] M. H. Azmi and J. Yuan, "Design of multi-edge type bilayer-expurgated 
LDPC codes," IEEE International Symposium on Information Theory 
(ISIT), pp. 1988-1992, Jun. 2009. 

T. M. Cover and A. A. E. Gamal, "Capacity theorems for the relay 
channel," IEEE Trans. Inf. Theory, vol. 25, no. 5, pp. 572-584, Sep. 
1979. 

T. J. Richardson and R. L. Urbanke, "Multi-edge type LDPC codes," 
Available Online: http://lthcwww.epfl.ch/papers/mutiledge.ps 
A. Chakrabarti, A. de Baynast, A. Sabharwal, and B. Aazhang "Low 
density parity check codes for the relay channel," IEEE J. of Set Areas 
Commun., vol. 25, pp. 280-291, Feb. 2007. 

C. Li, G. Yue, M. A. Khojastepour, X. Wang, and M. Madihian, "LDPC- 
coded cooperative relay systems: performance analysis and code design," 
IEEE Trans. Commun., vol. 56, no. 3, pp. 485-496, Mar. 2008. 
C. Li, G. Yue, M. A. Khojastepour, X. Wang, and M. Madihian, "LDPC 
code design for half-duplex cooperative relay," IEEE Trans. Wireless 
Commun., vol. 7, no. 11, pp. 4558-4567, Nov. 2008. 
R. Ahlswede, N. Cai, S.-Y. R. Li and R. W. Yeung, "Network infomation 
flow," IEEE Trans. Inf. Theory, vol. 46, no. 4, pp. 1204-1216, Jul. 2000. 
X. Bao, J. Li, "Adaptive network coded cooperation (ANCC) for wireless 
relay networks: matching code-on-graph with network-on-graph," IEEE 
Trans. Wireless Commun., vol. 7, no. 2, pp. 574-583, Feb. 2008. 
J. Kim, S. Park, J. Kim, Y. Kim, H. Song, "Joint LDPC Codes for Multi- 
User Relay Channel," Fourth Workshop on Network Coding, Theory and 
Applications (NetCod), pp. 1-6, Jan. 2008. 

S. Brink, "Convergence behavior of iteratively decoded parallel concate- 
nated codes," IEEE Trans. Commun., vol. 49, no. 10, pp. 1727-1737, Oct. 
2001. 

S. Brink, G. Kramer, and A.Ashikhmin, "Design of Low-density parity- 
check codes for modulation and detection," IEEE Trans. Commun., 
vol. 52, no. 4, pp. 670-678, Apr. 2004. 

G. Liva and M. Chiani, "Protograph LDPC codes design based on EXIT 
analysis," IEEE Global Telecommunications Conference (GLOBECOM), 
pp. 3250-3254, Nov. 2007. 

"A fast and accurate degree distribution optimizer for LDPC code 
ensembles," Available Online: http://lthcwww.epfl.ch/research/ldpcopt 
S. R. Kollu and H. Jafarkhai, "On the EXIT chart analysis of low- 
density parity-check codes," IEEE Global Telecommunications Confer- 
ence (GLOBECOM), pp. 1131-1136, Nov. 2005. 



[2] 



[4] 

[5] 
[6] 

[7] 

[8] 

[9] 
[10] 

[11] 

[12] 

[13] 

[14] 

[15] 
[16] 



